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Abstract. The CELESTE experiment uses the heliostats of an old solar farm in 
the French Pyrenees to detect gamma ray air showers by the atmospheric Cerenkov 
technique. CELESTE has been operating since November 1999 with an array of 40 
heliostats fully instrumented with 1GHz flash ADCs. Significant advances have been 
made in the detector simulations and in the data analysis techniques. We report here 
on results from recent observations of the Crab nebula above an energy threshold of 
50GeV. The results and simulations illustrate the current sensitivity of the experiment. 



INTRODUCTION 

CELESTE, described in detail in the experimental proposal [1] and in [2-4], 
uses the mirrors of a former solar electrical plant at the Themis site in the French 
Pyrenees. We have been operating since October 1999 with forty, 54 m 2 heliostats 
which direct the Cerenkov light, via secondary optics at the top of a 100 m tall 
tower, towards PMTs instrumented with flash ADCs. The "standard candle" of 
TeV astronomy, the Crab nebula, has been our primary target during this initial 
phase. 



SIMULATION RESULTS 

We have developed a simulation of CELESTE for use with standard air shower 
generators. During the 1999-2000 Crab observing season we collected 13 hours 
of data with both CAT and CELESTE which can be used to partially check our 
simulations. 20% of the CELESTE events, corresponding to ~30% of CAT events 
are common, as identified by their arrival time measured with GPS clocks at both 
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FIGURE 1. Simulation results. Left: The impact parameter (R) distribution as measured by 
CAT for real and simulated gamma rays. Right: The simulated response of CELESTE to an E~ 2 
power law gamma ray spectrum for a source at the position of the Crab at transit. 



experiments. Analysis of the CAT dataset [5] results in an excess of 1268 gamma 
events over a background of 3131 hadrons. The same analysis applied only to the 
common events gives 418 gammas over 526 hadrons. Thus, a CELESTE trigger 
doubles the signal to noise ratio in the CAT data, although it does not improve the 
significance as the data sample is smaller. CAT measures the shower impact param- 
eter with a resolution of ~5 m. Fig. 1 shows this reconstructed impact parameter 
for simulated data, and for excess events from the common Crab dataset. The sim- 
ulations reproduce the data, showing that the effective surface area for CELESTE 
in the energy region of the common CAT-CELESTE events is understood. 

Also shown in Fig. 1 is the energy distribution of simulated events for an input 
E~ 2 differential gamma ray spectrum. At the raw trigger level the energy threshold, 
defined as the peak differential gamma ray flux, is ~30 GeV. We have found it 
necessary to apply a software trigger in order to remove the effects of night sky 
background differences between the ON source and OFF source control region, 
increasing our energy threshold to ~50 GeV. We also require > 10 FADC peaks with 
an amplitude > 25 digital counts to ensure the shower parameters are measurable. 

We have used the simulations to develop our analysis techniques. The FADC 
peaks are fitted with a function modelled on the single photo-electron pulse shape 
and the reconstructed charge and timing information is used to derive the pa- 
rameters of the shower. The homogeneity of the light distribution on the ground 
provides the best discrimination between hadronic and gamma ray showers. We de- 
fine a parameter <J grp , calculated by summing the FADC signals from 8 heliostats in 
5 groups, corresponding to the hardware analogue sums used to trigger the experi- 
ment, then calculating the variance of the amplitude of the 5 pulses and normalising 
this variance to the mean amplitude. The distribution of a grp for real data and for 
the simulations is shown in Fig. 2. A cut at a grp <0.25 is expected to reject 80% 
of the remaining hadronic background and retain 70% of the gammas. 

Also shown in Fig. 2 is the distribution of the angle, 9, between the shower axis 
and the direction to the source reconstructed using two points: the shower impact 
parameter at ground level, defined by the barycentre of the charge distribution over 
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FIGURE 2. Analysis parameters. Left: The distribution of a grp for simulated gammas, protons 
and real OFF source data. Right: The same for the shower axis angle, 9. 



TABLE 1. The number of events, excess, significance and signal to 
background ratio (S/B) after each cut for the Crab dataest. 





ON 


OFF 


Excess 


Significance 


S/B 


Raw 


850 922 


845 021 


5 901 


4.5 


0.7% 


Software Trigger 


429 533 


424 881 


4 652 


5.0 


1.1% 


A poaks > 10 


402 452 


398 242 


4 210 


4.7 


1.1% 


Cgrp < 0.25 


84 499 


81937 


2 562 


6.2 


3.1% 


9 < 7 mrad 


49 244 


47466 


1778 


5.7 


3.7% 



the heliostats, and the position of the shower core relative to the pointing direction 
at the mean altitude of maximum emission for gamma showers in the CELESTE 
energy range (11 km above the altitude of the experiment). The latter point is 
calculated by fitting a sphere to the measured arrival time of the Cerenkov wavefront 
at each heliostat. The small field of view of CELESTE (10 mrad fullwidth) reduces 
the observable differences between hadron and gamma showers; however, a cut 
at 9 <7 mrad is expected to reject 20% of the remaining protons for 10% of the 
gammas. The effect of the analysis cuts on the energy distribution of the gamma 
showers is shown in Fig. 1. 

CELESTE OBSERVATIONS OF THE CRAB NEBULA 

The Crab Nebula, as a strong steady source of VHE gamma rays [7] , is the obvious 
first target for CELESTE. We have collected 11. 5h of ON source and corresponding 
OFF source data with all heliostats tracking a point 11 km above the altitude of 
the experiment. Fig. 3 shows the distribution of 9 after all cuts, with an excess of 
events at small 9 as expected. Table 1 shows the events remaining at each stage of 
the analysis. The final significance of the result is 5.7 a 

We have calculated the normalisation factor to an assumed E~ 2 differential spec- 
trum from the measured gamma ray rate, giving d§/dE = 

(2.5 ± 0.5 stat ± 1.2 syst ) x 10" 8 /E 2 photons cm" 2 GeV" 1 s" 1 




FIGURE 3. Results on the Crab. Left: The difference in the distribution of 9 between ON and 
OFF source observations. Right: The Crab spectrum, including the preliminary CELESTE result 
and a model for the synchrotron and inverse Compton emission (after [6]). 



The CELESTE measurement is shown in Fig. 3 along with the results from other 
experiments. This is a preliminary estimate, as we are still quantifying systematic 
errors in various areas of the analysis, particularly in the comparison between the 
predicted and observed efficiency of analysis cuts. 

We have also made some observations of the Crab using different tracking meth- 
ods. A short (2h 10m) exposure using a double pointing method, wherein half of 
the heliostats point at 11 km altitude and the other half point at 25 km, appears 
to double our sensitivity by allowing us to sample a greater part of the shower. 



CONCLUSIONS 

The CELESTE experiment has detected the Crab nebula with high significance 
(5.7 a). Simulations have been used to develop our analysis procedure and also 
to estimate our energy threshold to be ~ 30 GeV at the raw trigger level and 
~ 50 GeV after the analysis cuts. CELESTE will be upgrading to 54 heliostats 
shortly. This, in combination with improvements to the electronics and the double 
pointing strategy should greatly improve our sensitivity for the next Crab season. 
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